A Chiral Paramagnetic Skyrmion-like Phase in MnSi 
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We present a comprehensive study of chiral fluctuations in the reference helimagnet MnSi by 
polarized neutron scattering and Neutron Spin Echo spectroscopy, which reveals the existence of 
a completely left-handed and dynamically disordered phase. This phase may be identified as a 
spontaneous skyrmion phase; it appears in a limited temperature range just above the helical 
transition Tc and coexists with the helical phase at Tc- 



Chirality is ubiquitous in nature and of fundamental 
importance both on the microscopic level and in our ev- 
eryday hfe. The break of symmetry between right and 
left manifests itself in parity violation, governs biolog- 
ical structures such as DNA and can also be experi- 
enced in the organisation of our own body. In mag- 
netism, chirality is evident in solitons [1], systems with 
geometric frustration [2] and metallic systems with non- 
centro-symmetric lattice structures, where the resulting 
anti-symmetric Dzyaloshinski-Moriya (DM) interactions 
[3, 4] introduce a parity breaking term in the Hamilto- 
nian [5]. The DM term has the form Af x (V x M) and 
is more than a perturbation giving rise to the peculiar 
canted magnetic arrangements found in high tempera- 
ture superconductors [6] or the cycloid spin structures 
in multiferroics [7, 8]. In the non-centrosymmetric weak 
itinerant-electron ferromagnet MnSi, DM induced chiral- 
ity comes in close interplay with Fermi liquid behavior 
and quantum fluctuations [9] . The Hamiltonian of MnSi 
comprises three hierarchically ordered magnetic interac- 
tion terms with well separated energy scales [10], which 
allow to distinguish between different contributions. The 
strongest ferromagnetic exchange interaction aligns the 
spins, the weaker chiral Dzyaloshinski-Moriya (DM) term 
twists them into a helix and the weakest Anisotropic Ex- 
change (AE) or crystal field term pins the helix propaga- 
tion vector f along the (111) crystallographic directions. 
The helical order appears below Tc w 29 K. It is a left- 
handed helix with a period of £ ~ 175 A (r w 0.036A~^) 
and all magnetic moments perpendicular to the helix vec- 
tor [11]. 

In this letter we concentrate on the chiral correlated 
paramagnetic or spin liquid phase of MnSi just above Tc, 
where intense diffuse neutron scattering spreads homoge- 
neously over the surface of a sphere with radius t. This 
unusual feature emerges as a ring on the two-dimensional 
small angle neutron scattering patterns and the rings re- 
duce to half-moons if the beam is polarized. This is illus- 
trated by figure 1, which reproduces spectra from [12]. 
Numerous theoretical studies were devoted to explain 
this phase invoking possibilities such as unpinned helical 



order [12, 13] or condensation of chiral order parameters 
[14]. Recent local mean-field calculations assuming the 
hierarchical hamiltonian of MnSi show that the helical 
phase is preceded by a disordered phase with skyrmion- 
like short range order similar to the partial order in liquid 
crystals [15], which sets in at Tc ~ Tc + IK (see supple- 
mentary information of [15]). Skyrmions are solutions of 
the non-linear field theory introduced by Skyrme [16] to 
bridge the gap between waves and particles in the wave- 
particle duality concept. Skyrmionic, vortex-like, states 
and lattices were first predicted in non-centrosymmetric 
magnetic systems by A. N. Bogdanov [17]. In MnSi a 
skyrmion lattice was recently identified under magnetic 
field [18]. At zero magnetic field skyrmions would result 
from a vortex-like short range order of the isotropically 
twisted spin liquid phase (see insert on Figure 2b). 

Spin correlations and fiuctuations of MnSi have been 
extensively studied by triple axis spectroscopy (TAS) at 
Brookhaven in the mid-80s [19] and more recently with 
polarized neutrons [20]. These studies, however, did not 
have the resolution needed to see fiuctuations close to the 
helical transition. We opted for neutron spin echo spec- 
troscopy (NSE) [21], which reaches the highest energy 
resolution (<1 neV) among all neutron scattering tech- 
niques. We also combined NSE with spherical neutron 
polarimetry (SNP) [22] to the polarimetric NSE tech- 
nique [23] for a rigorous analysis of chiral fiuctuations. 
The starting point of our study is the formalism devel- 
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FIG. 1: Polarized neutron small angle scattering patterns 
measured by [12] (a) slightly below (Tc-0.1 K) and (b) above 
(rc-|-0.2 K) the helimagnetic transition of MnSi. The vectors 
Pa mark the polarization of the incoming neutron beam. 
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oped almost simultaneously by Blume [24] and Maleyev 
[25] in the early '60s. The magnetic scattering cross- 
section for thermal neutrons may be written as: 



M± ■ Ml + P' ■ 9(Mx X M^) 



(1) 



where P' is the polarization vector of the incident beam 
and the magnetic interaction vector M± is the projection 
of the magnetic structure factor onto a plane perpendic- 
ular to the scattering vector Q. The (imaginary) vec- 
tor product is the chiral term characteristic of the spiral 
structure. 

In MnSi the helical arrangent of the magnetic moments 
is described by two orthogonal vectors of the same am- 
plitude, which are perpendicular to the helix vector r. 
At the helical Bragg peaks fm the sample orientation is 
such that for a left handed helix the cross-section is max- 
imum for P' = Q and completely extinct for P' = —Q 
[24]. This extinction is responsible for the characteris- 
tic half-moon small angle scattering patterns seen with 
polarized neutrons above Tc and gives a very precise 
determination of the weight of the dominant chiral do- 
main 7] = |M_L X M^|/(Mx • Ml). T] can also be de- 
termined from the polarization vector of the scattered 
neutron beam P: 



Pa = -P'{AJ^- Ml) + 2di{{P' ■ Ml)Mj_). 
... -3(Mx X MX) 



(2) 



The first two terms lead to the conventional tt rotation 
of the polarization vector around M± . The chiral term 
creates a polarization (anti)parallel to Q independently 
from the incoming beam polarization. Complete knowl- 
edge of Eq. 2 implies complete knowledge and control of 
P' and P independently from each other. This is achieved 
with spherical neutron polarimetry using Cryopad, which 
brings P' and P either along x = Q, z (perpendicular to 
the scattering plane) or y, which completes the right- 
handed cartesian set [26] . This procedure determines the 
polarization matrix: 
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with {a, (3) e {x,y,z} (3) 



where P is the polarization tensor of the non-chiral terms 
and P^ the polarization created by the chiral term. At 
fill, X II T and only the first row of the matrix is non- 
zero. The ideal matrix is given in Table I, where C de- 
termines the chirality of the helix: C = -|-1 for right and 
C = — 1 for left handed chirality. rj measures the fraction 
of the dominant chiral domain, = 1 for a single domain 
and = for a disordered state or equally populated chi- 
ral domains. 

The polarimetric measurements were done with Cry- 
opad on the NSE spectrometer IN15 of the ILL. With 
the polarimetric NSE configuration it was possible to 



TABLE I: Ideal and measured polarization matrix Pa,/3 at 
fill. For details see text. 
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FIG. 2: Intensity at the position of the fm helical peak (a), 
chiral fraction of the dominant left handed domain (b) and 
characteristic time to of the fluctuations at fm and fno (c). 
In (b) the red circles are deduced from the extinction (Px,x 
eq. 1), the blue squares from the non-diagonal terms Vy,x 
and Pj,:c and the green dotted line is calculated from ^m 
and eq. 4. The error bars are calculated from the counting 
rates and do not include systematic errors. For the sake of 
comparison between difl^erent cryogenic setups the data are 
plotted against T — Tc. The lines are guides for the eyes. 



determine the matrix Va,f3 and at the same time ac- 
cess the intermediate scattering function I(Q,t) of each 
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non-zero term. These measurements were complemented 
by "conventional" NSE measurements on INll (ILL). The 
measurements were done at A=9 A (IN15) and 5.5 A 
(INll) with AA/A=15% FWHM and a Q-resolution of 
~ 0.005 A"'^ FWHM. Most of the measurements were 
carried out at fm- Additional spectra were recorded on 
INll around tuq (the position of tuq is shown on the 
SANS patterns of figure 1). It was also possible to ex- 
tract precise information on the q = Q — t dependence 
of the relaxation by analyzing slices of the detector at 
constant distance from the 000 point. The MnSi single 
crystal was the same as in previous experiments [12] with 
the B20 lattice structure {P 2i 3) and a lattice constant 
of 4.558 A. 

The transition temperature Tc = 29.05 ± 0.05 K is 
defined from the sudden increase by more than one order 
of magnitude of the neutron intensity at rm (figure 2a), 
which confirms the weak first order transition with strong 
second-order features seen on specific heat [27]. 

The Ornstein-Zernike analysis of the q dependence of 
the intensity around -fm and fno leads to correlation 
lengths 4iii and ^no, which are very close to previous re- 
sults [12] and are shown on the insert of figure 2a. These 
quantify the extent of correlations along the helix vector 
f. 4iio stays significantly smaller than the pitch of the 
helix and shows a maximum slightly above Tc, whereas 
^111 diverges and reaches the pitch of the helix (dashed 
line) at ~ Tc-|-0.05 K triggering the first order transition. 

The polarization matrix was measured with an inci- 
dent beam polarization of 97% and the data were cor- 
rected for the background signal. At fm the calculated 
and measured matrices of the helical phase at 25 K are 
in excellent agreement with each other (Table I) confirm- 
ing the stabilisation of a single-domain left-handed helix 
with r] = l,C = ~1- When increasing the temperature, 
the polarization matrix Pq,/3 is unaffected by the first or- 
der transition and remains unchanged up to Tc Tc + l 
K. Also the extinction of the scattering when P' = —t 
(eq. 1) is complete up to Tc- All values of t], plotted in 
figure 2b, reveal the existence of a new phase between Tc 
and Tc ■ This is a completely single domain phase even 
though the correlations are significantly shorter than t. 
at Tc the correlation length does not exceed ~£/4. This 
result contrasts with the common experimental [28] and 
theoretical [12] picture based on the hierarchy of the in- 
teractions, according to which MnSi behaves as a ferro- 
magnet at distances shorter than £ and for ^ < ^. In fact 
our experimental results show a more complex situation 
and the slow decrease of the chiral fraction 1] displayed 
in figure 2c spreads over a very wide temperature range. 
Even at 50 K (Tc +20 K) where correlations do not ex- 
ceed 7 A [19], we found r] ~ 0.203 ± 0.024. 

A mean field description of paramagnetic fluctuations 
in a ferromagnet in the presence of DM interactions was 
given in [12] and renders the main features of MnSi: ex- 
tinction of the scattering at P' = — f or intensity maxi- 



mum at the sphere Q = t (flgure 1). The theory also 
gives a simple relation between rj and the correlation 
length 4: 

ri = 2QT/{Q^+T^ + l/e) (4) 

From the measured ^m values we calculated the dotted 
green line on figure 2, which is significantly lower than the 
experimental data. Consequently, the simple assumption 
of an unpinned or fiuctuating helical phase is not suffi- 
cient to describe the experimental result. The high values 
of T] refiect a partial or short range order underlining the 
non-trivial character of the spin liquid phase of MnSi, in 
particular of the phase between Tc and Tc ■ 

In the following we will discuss the results of our 
NSE measurements. NSE spectroscopy gives the nor- 
malized intermediate scattering function I{Q, t) = 
S{Q,t)/S{Q,0), which is the time-energy Fourier trans- 
form of the scattering function S(Q,w). The fiuctuations 
of the single domain chiral spin liquid phase were mea- 
sured with standard NSE. The fluctuations for P' _L rm 
(Pz,x or Py^x) and P'llrm (Px,x) at the fm position were 
measured with polarimetric NSE. Figure 3 shows the in- 
termediate scattering function for the non-diagonal term 
IPz,x- We flnd exponential relaxations superimposed on 
an elastic contribution due to the Bragg peak, which were 
fltted by the function I{Q,t) = aexp{—t/to) -|- (1 — a). 
The elastic fraction (1 — a) evolves from ~ 20% to 100% 
within 0.2 K following the fast increase of the intensity in 
flgure 2a. The onset of the Bragg peak masks the decay 
of I{Q,t) at Tc- On the other hand, at fno there are 
no Bragg peaks and fluctuations may be seen also be- 
low Tc ■ The measurements were extended up to Tc + 2 
K, where they agree with the previous TAS [19, 20] and 
preliminary resonance NSE results [29]. 

All values of to are plotted against temperature in flg- 
ure 2c. The existence of flnite relaxation times at these 
temperatures rules out the interpretation in terms of 
quasi-static correlations suggested by [9]. The measured 
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FIG. 3; Dynamic correlations at the position of one of the 
helical Bragg peaks fm measured by polarimetric NSE for 
the ¥^^x term of the polarization matrix (Table I). 
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relaxation times are rather typical for correlated systems 
such as magnetic nanoparticles [30] or ferromagnets at 
Tc [31]. All t() (at Till and tho) fall on top of each other 
within experimental error bars: the direction of the in- 
cident polarization vector or the particular position on 
the sphere with radius r have no effect on the observed 
dynamics. The fluctuations are governed by the corre- 
lation length at rno without any sign of slowing down 
at Till, i-e. at the position where the Bragg peak ap- 
pears. Consequently, the fluctuating chiral spin liquid 
phase does not develop to the helical phase. The two 
phases coexist close to Tc and may be decomposed 
as: 6ii = a^fiuct + (!-«) Astatic with ^fi^ct = 6io re- 
flecting the fluctuating spin liquid fraction a and Astatic 
the static component of the helical phase. It has been 
suggested that fluctuations may drive the transition in 
MnSi to first order [10, 32, 33]. Our present high reso- 
lution data rule out any interplay between fluctuations 
and the first order helical transition. 

Dynamic scaling relates the characteristic relaxation 
time to the correlation length: to cc , with z the dy- 
namic exponent. Below Tc we found to oc ^fiuct leading 
to 2; « 1, an unconventional value much lower than in 
3D ferromagnets (2=5/2) or antiferromagnets {z=3/2). 
Such a low value could refiect a locally reduced dimen- 
sionality with the spins pinned at the local helical planes. 

Above Tc the Bragg peaks disappear and the cor- 
relation length becomes considerably shorter than the 
pitch of the helix. In spite of that all spins remain 
pinned into the local helical plane up to Tc'- This is 
exactly the theoretically predicted temperature of forma- 
tion of a skyrmionic ground state in MnSi [15]. It is also 
the temperature range of a broad maximum in the spe- 
cific heat, which refiects strong short range correlations 
[27]. Besides specific heat, thermal expansion, sound ve- 
locity, sound absorption as well as resistivity measure- 
ments point towards a non trivial spin liquid phase [34] . 
Skyrmions account for the inherent topology of this sin- 
gle domain spin liquid phase, which gives a scattering 
with rotational symmetry. Skyrmionic textures are ener- 
getically favourable at short distances, which privileges 
structures with magnetic moments decreasing with in- 
creasing distance from the vortex center [15]. In the spin 
liquid phase the mean magnetic moment decreases con- 
tinuously with increasing distance d following the typ- 
ical correlation function ~ exp(— d/^)/d. The correla- 
tion length ^fiuct = Clio would, therefore, give a natural 
measurement of the skyrmion radius, which would not 
exceed £/2 at Tc stabilizing the single domain spin liq- 
uid phase. That skyrmionic solutions are energetically 
more favourable at short lengths would explain the slow 
decrease of rj above Tc and its finite value even at 50 K. 

In summary, these first polarimetric NSE measure- 
ments provide the missing link in getting a consistent pic- 
ture of the spin liquid phase of MnSi above Tc- High res- 
olution neutron spin echo spectroscopy and polarimetry 



evidence the existence of a disordered skyrmion phase, a 
completely chiral single domain and strongly fiuctuating 
new state of matter in a very narrow temperature range 
of ~ 1 K above the helical phase of MnSi. 
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